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We investigate the relative phase between two weakly interacting ID condensates of bosonic atoms 
after suddenly switching on the tunnel-coupling. The following phase dynamics is governed by the 
quantum sine-Gordon equation. In the semiclassical limit of weak interactions, we observe the 
parametric amplification of quantum fluctuations leading to the formation of breathers with a finite 
lifetime. The typical lifetime and density of the these 'quasibreathers' are derived employing exact 
solutions of the classical sine-Gordon equation. Both depend on the initial relative phase between 
the condensates, which is considered as a tunable parameter. 



Dynamical instabilities can amplify spatial field fluc- 
tuations drastically. If the instability provides sufficient 
energy, even quantum zero-point fluctuations can trig- 
ger the formation of macroscopic field patterns. For in- 
stance, in some cosmological scenarios of the inflationary 
stage of the universe, following a slow-roll, a scalar in- 
flaton field performs oscillations around the minimum of 
the corresponding inflaton potential. Thereby, resonant 
spatial fluctuations are amplified parametrically. For a 
rather generic class of potentials these can end in long- 
lived, local concentrations of energy, so-called 'oscillons' 
(e.g.,[l-4j). 

In this work, we argue that analogous non-equilibrium 
phenomena should be observable in experiments with 
a pair of weakly interacting quasi- ID clouds of cold, 
bosonic atoms [5j|6]. After suddenly turning on the tun- 
nel coupling between the condensates, the dynamics of 
the relative phase field <\> is governed by the integrable 
quantum sine-Gordon model |7J 
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where the 'mass' m depends on the tunnel amplitude, 
such that m(t) = mQ(t) for this quench. The phase has 
been rescaled, and for weak interactions j3 <C 1. 

The sine-Gordon model (SGM) is one of the 
most prominent prototypical models of low-dimensional 
condensed-matter systems. Currently, quenches in the 
SGM are under intense investigation (e.g., [8UTT]). For 
instance, only recently the amplification of density inho- 
mogeneities after a quench to the special Luther-Emmery 
point (f3 2 = 47r), was demonstrated |12| . 

Here, we propose that the spatially averaged value <E> 
of the relative phase field is tuned to some value <3>o right 
before the quench (e.g., this might be achieved by slightly 
tilting the transversal double- well potential confining the 
BECs, cf. Fig.[TJi). The subsequent phase dynamics can 
be directly observed in matter wave interference experi- 
ments [5j[6]. 

At short times, will perform Josephson oscillations 
(Fig.JTJ)) according to d 2 &/dt 2 = -j3~ lr m 2 sin/3<£. How- 




Figure 1. a) Proposed experimental protocol. A single quasi- 
1D condensate is split and the global phase is tuned to $o- 
Switching on the tunnel-coupling, the dynamics of the relative 
phase obeys the quantum sine-Gordon equation, b) At short 
times, the global phase performs Josephson oscillations c) Sin- 
gle run of TWA (see main text) with f3&o — 0.257T. Spatial 
quantum fluctuations are amplified parametrically. Eventu- 
ally, the non-linearity of the sine-Gordon equation kicks in 
and 'quasibreathers' (breathers with a finite lifetime) form. 



ever, these are linearly unstable [T3| fl4 ] in the presence of 
inhomogeneous quantum fluctuations which are paramet- 
rically amplified at certain wave lengths (Fig.[I]:). As we 
will show in this paper, due to this modulation instability 
at some later point the dynamics becomes fully nonlinear 
and one observes the formation of sharply localized and 
oscillating patterns in the phase-field (Fig. [I]:). It will 
be shown below that they can be related to particular 
exact solutions of the classical SGM obtained from the 
Backlund transformation (following |15|). These 'qua- 
sibreathers' (QBs) have a finite lifetime (in contrast to 
the well-known breather solutions), and we find good ev- 
idence that a quasi-equilibrium steady-state with a finite 
density of such excitations develops at long times. 
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For the considered quench, due to the intrinsic insta- 
bility, a simple approximation of the SGM in terms of 
non-interacting, massive phonons is invalid, even though 
f3 <C 1. Instead, we employ the truncated Wigner- 
approximation (e.g., [16J) (TWA). The basic idea be- 
hind TWA is to simulate classical field equations, but 
with quantum-mechanical fluctuations as stochastic ini- 
tial conditions. In the limit of /? —> 0, TWA is expected to 
become reliable as it can be shown that (3 plays the role of 
an effective Planck's constant [16J. It describes correctly 
the linear dynamics during the parametric amplification. 
Once QBs form, occupation numbers are already large 
such that a semiclassical description (provided by TWA) 
should continue to remain valid. The great advantage of 
TWA is that it serves snapshots (Fig. [I]:) of the phase 
field <j)(X) t) run-by-run, which can be compared directly 
to the experimental observations. 

In the following, we introduce the model and demon- 
strate that the proposed setup should be well within the 
reach of present experiments. After a discussion of the 
numerical findings obtained from TWA, we will intro- 
duce the analytical 'quasibreather' solutions of the clas- 
sical SGM. We demonstrate that these solutions are well 
suited to explain the main physical features as predicted 
by TWA. Eventually, we argue that a statistical analysis 
of </>(x,t), experimentally obtained at a single time t per 
run, could reveal the distinctive signature of QBs. 

Model - It was shown [7] that on scales larger than the 
condensate healing length the dynamics of the relative 
phase between condensates <p = (0i — (/>2)/V^ (and of 
the density field ft, fulfilling [Il(x),(j>(x')] = ihS(x - x')) 
is governed by the quantum SGM. In the following, we 
use the notation in |7J. After rescaling the fields ((J) \-> 
^7r/K^) and n \-> y/K/irfl), the Hamiltonian reads 



H 



dx 



{d x <t>y 



P 2 



dx cos (3<j). (2) 



The repulsive short-range interaction, characterized by 
the Luttinger parameter K, enters (3 = y27r / K. For 
repulsive bosons K G [l,oo]. In the limit of small 
interactions [7J (K — >► oo), K is connected to the Lieb- 
Liniger parameter 7 = rriBg/h 2 po via K = tt/^/j (where 
g is the interaction strength, m# the mass of the bosons 
and po the mean density). Furthermore, there are the 
sound velocity v s = vp/K (with vp = Httpo/tub), and 
the effective mass m = ^2t±po/3, where the tunnel- 
amplitude t± between the condensates enters. The rele- 
vant length scale in the sine-Gordon model, determining 
the width of breathers and solitons, is set by ^hv s /m. 
The sine-Gordon model should be a sound description of 
the system as long as y/v s H/m ^> ~ 1/po-^/j. Finally, 
the Josephson frequency is given by ujj = m^/v s /h [T7] , 
determining the period of the Josephson oscillations for 
[3&o < 1. In the following, we set v s = h = 1 . 

Experimental Realizability. - In the limit of weak inter- 
actions <C 1), the stiffness of the condensate becomes 




Figure 2. a) Space-time plot of the phase- field for a single 
simulation run. b) Spatial cut through the quasibreather 
highlighted in a) over one breather period (thin blue line; 
thick blue line: filtering out irrelevant short wavelength fluc- 
tuations), c) Temporal cut through the emerging QB taken 
at its center. 5<j) = <j) — $ is obtained by subtracting the 
spatially averaged field $ = \ J Q L dx cj). One observes the ex- 
ponential amplification of initial quantum fluctuations. Here, 
f3$ = 0.37T. 



small. To avoid the breakdown of the hydrodynami- 
cal description, the mean condensate density po must 
be large enough. In fact, the condition ^v s h/m ^> £^ 
translates into a lower bound for f3 (and therefore into 
a lower bound for the interaction strength): j3 2 ^> 
4po 1 VtjjriB/h. Within the present experimental setups 
the condensate density can be widely tuned. E.g., in [T5| 
the ID density of Rb r ~ 
ranges from po = 3 
frequency uj± = 2ir x 4kHz, with g = 2hw±a s and 
a s = 5.31 • 10 _9 m, we have K = 9 — 52 and therefore 
small values of the sine-Gordon parameter f3 = O(10~ 1 ) 
are well achievable. Eventually, the tunnel-amplitude can 
be tuned between t±/h = [O(10 2 ) - O(10 4 )]/s [5J. Thus, 
for large densities po = 100/im -1 the lower bound on f3 
is O(10 _1 ) as well. In all numerical simulations, we use 
j3 = 0.1. We conclude that after some fine-tuning of the 
experimental parameters, the proposed setup should be 
within reach. 



87 atoms (with m B = 1.44- 10~ 25 kg) 
100/im -1 . For a transverse trap 



Numerical results - Before the quench, the phase field 
(j) is given by the offset $0 plus the small zero-point fluc- 
tuations of the free theory [Eq. ([2| with m = 0]. Within 
TWA, the fluctuations of the (j) and n modes are ini- 
tialized according to their Gaussian Wigner-distribution 
[16], omitting the q = mode. 

This initial state should be experimentally achievable. 
One should start from the ground state at strong tunnel 
coupling (m ^> 1), where the relative phase <j>(x) = 0. 
During a slow linear rampdown of the tunnel-amplitude, 
all modes (j) q ^o will follow their time-dependent ground- 
state. It can be shown that the global phase $ at the end 
of this process obeys f3 2 (& 2 } ~ l/£hPo, which remains 
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small as long as the number of bosons within the healing 
length is large. An additional potential tilt will produce a 
fixed phase offset <£> ~ <3>o- We will focus on rather small 
f3<f>o < 1. For — > 7r, in addition to breathers one 
observes the formation of solitons. Finally, we note that 
this 'phase-tuning' becomes impossible for infinite sys- 
tems. This is due to the logarithmic divergence of phase 
fluctuations with system size L in the massless ground- 
state, i.e., P 2 (5(j) 2 } ~ p 2 \nL. However, these are largely 
suppressed in the semiclassical limit ft <Cl. We always 
choose the system size such that the initial overall rela- 
tive phase is well denned. 

In Fig.[2j the phase field for a single run of TWA 
is shown. At small times, the field is dominated by 
$>(t) performing ordinary Josephson oscillations. At this 
point, the inhomogeneous part 8(j)(x,t) of (j) = $ + 5<fi 
can be treated as a small perturbation. The modes obey 
d 2 (j) k + (k 2 +m 2 cos[/3$(£)])0 fe ~ for all k ^ 0, i.e., these 
are phonons with a periodically modulated mass. Modes 
with \k\ G [0, m sin |/J3>o/2 |] are parametrically amplified 
yielding 4>k(t) ~ e Tkt . The amplification rates T k for this 
linear regime (here displayed for /3<l>o < 1) (see Fig.[3^i) 

2T k ~ ^sin 2 (^ /2)-P/m 2 , (3) 

were found in [13J, neglecting the damping of the driving 
<3>-mode. After some time, nonlinear interactions between 
the amplified modes become important and lead to the 
formation of sharply localized oscillating structures (here 
denoted as 'quasibreathers'), which constitute the main 
phenomenon discussed in our paper. Once these localized 
oscillations get out of phase with respect to the back- 
ground oscillations of <£, their energy is depleted again. 
One arrives at a steady-state, where QBs are randomly 
created and decay, with a typical lifetime f and a mean 
spatial distance A (cf. Fig.[2J). All these statistical quanti- 
ties crucially depend on the initial value <3>o. It turns out 
that these localized modulations in the stochastic phase 
field can be connected to certain exact solutions of the 
SGM. These solutions, to be discussed in the following, 
are standing breathers riding on top of a homogeneous 
and oscillating background. We will demonstrate that 
this set of solutions is well suited to describe the numer- 
ical observations and provides analytical insight into the 
dependence of A and f on $o- 

Quasibreathers. - A Backlund transformation (see for 
instance [15J) allows to 'add' (ant i) solitons to a given 
solution of the SGM, here taken to be the spatially ho- 
mogeneous solution $>(t) with $(0) = &o an d $(0) = 0. 
Adapting the approach of [15J, we obtain quasibreather- 
solutions by adding a (standing) soliton and the corre- 
sponding antisoliton to A standing SGM-breather is 
characterized by its amplitude (e.g., [19J). The solutions 
here depend in addition on the amplitude of the under- 
lying background oscillations <3>o (furthermore, there is 
a minor dependence on the precise value of the relative 
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Figure 3. Phononic zero-point fluctuations are amplified para- 
metrically driven by the overall phase $ performing Joseph- 
son oscillations, a) Corresponding amplification rate Tk of 
phonons with wavenumber k for /3&o/tv G 0.05 — 0.4 accord- 
ing to Eq. Q (dashed line). Solid lines show the rate obtained 
from the 'multibreather' solution Eq. j5|. b) Plot of a (spa- 
tially periodic) 'multibreather' solution reducing to a driven 
phonon with wavenumber k — mf3^o/^/8 for t/r% — > — oo. 
Cut A shows the parametric amplification of this phonon (B) 
driven by the oscillating background (we plot Sxk, substract- 
ing this background). They provide 'seeds' for the formation 
of Quasibreathers (C) with a lifetime r%. Here, f3&o = 0.3tt. 
c) Plot of a 'single quasibreather' (periodic in time) with 
/3<J> — 0.37T and ip = 1.02. d) Mean frequency Q = m — 
of a 'single quasibreather' (| smcp\ < | cos ^|^-|) depending on 
the amplitude of the background oscillations and the param- 
eter ip. For $o=0, one restores the unperturbed breather 
frequency Q=msm(p. 

phase Ao between breather and background at t = 0). 
They have the form 

4 

X(x, t) = p arctan t; $ , A , cp)] + (4) 

For $ = 0? the parameter (p G [0,7r/2] determines the 
unperturbed breather frequency rasing and max% = 
P~ 1 (2tt — 4(p). An explicit expression for Q and a de- 
tailed discussion of Eq. Q are given in the supplement. 
Here, we focus on the most relevant features in the limit 
0*0 < 1. 

The central observation is that by placing a breather 
on top of the oscillating background, its amplitude be- 
comes time-dependent (see Figs.[3]3,c). The background 
<fr(t) amplifies the breather while the frequency of the 
latter decreases. This is due to the fact that the effec- 
tive curvature of the cosine-potential decreases for larger 
field amplitudes. Eventually, both run out of phase and, 
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subsequently, the breather gets damped, 
phase-drift occurs at a frequency 



The relative 



ra^|cos 2 (/3$ /2) 



sm (f cos (/? 



1 + sin ip 



(5) 



As long as sin(/? < cos(/3$o/2), x describes a single 
breather whose amplitude is modulated with a period 
2ir/5Q (Fig. [3]:). This constitutes a stable QB solution of 
the SGM. As the background oscillates at frequency ra, 
one can understand Q = m — SQ as the mean frequency 
of the quasibreather (Fig.[3]i). 

For sin(/? — »> cos(/3<I>o/2), however, the period diverges. 
In fact, it turns out that for sirup > cos(/3<I>o/2), x 
is periodic in space rather than in time (Fig. [3)3). It 
describes a set of quasibreathers at distance 2ir/k = 



2nm 



cos' 



C8*o/2) 



■ sin 2 ip I 



-1/2 



These 'multibreather' solutions have remarkable prop- 
erties. For t — > — 00, they reduce to a phonon #x/e with 
wavenumber /c, superimposed on the oscillating back- 
ground. It is amplified parametrically, yielding 



Sxk oc e k C(t) coskx. 



(6) 



where ((t) is a periodic function and the amplification 
rate = 5Q((p(k),Q>o). This rate is in agreement with 
the results in Eq. ([3| for small /3$o- At later times, the 
characteristic breather peaks form (Fig. [3]:), exist during 
a time set by = 4/1^ and decay again for t — »• 00. Note 
that x therefore describes the prototypical formation of 
QBs out of fluctuations in a single mode 

Our numerical analysis shows (see below) that one can 
infer the properties of typical QBs observed within the 
stochastic TWA from these 'ideal' solutions x- In a given 
run, zero-point fluctuations in all modes are present. 
However, in the weakly interacting limit (/? <C 1), the 
parametric instability automatically filters out modes 
with k ~ fc, where k ~ m(3$o/y/8 denotes the maximally 
amplified mode. Thus, the typical distance between QBs 
is roughly given by A = 2ir/k. Although the strictly pe- 
riodic 'multibreather' solution x is not directly observed, 
individual QBs are well described by x as long as their 
amplitude is large enough compared to the noisy back- 
ground. Typical QBs decay after a lifetime f = (cf. 
Fig. [3^). When their amplitude is depleted down to the 
noise level, the exact solution ceases to be relevant. Then 
the background oscillations initiate the amplification pro- 
cess again, leading to the formation of new QBs. This 
explains the observed stochastic creation and annihila- 
tion of QBs at large times T^t ^> 1. 

Statistical Analysis. - These predictions agree well 
with the statistical analysis of the quench within TWA. 
For every run, we numerically track all QBs, showing 
their shape (/)qb(x) at the time of maximum amplitude 
(Fig.|4|i). A comparison of the mean QB-shape (4>qb) 
and x (with k = k) shows excellent agreement. Note 
that no fit parameter enters here. 
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Figure 4. a) Mean shape of QBs taken at their maximal 
value, obtained from TWA (blue dots) and full distribution 
of shapes (color plot). We track QBs numerically run- by- 
run in the interval t/f G [0.5,4] for P&o/tt = 0.3. Com- 
parison to the analytical solution x( x i 0) with k (blue, solid 
line; Ao = 0) shows almost perfect agreement. Inset) Var- 
ious /3$ /vr = 0.3,0.4,0.5 (from bottom to top). No fit- 
parameter enters, b) Correlation function C<^, averaged over 
time t/f e [2.25,3] (dots) with P$ /tt = 0.3,0.4,0.5 (from 
bottom to top). Locally, it agrees well with Eq. |7|. Inset 
shows 600(0, £), demonstrating that the field enters a steady 
state for t/f > 1. 



Finally, Fig.[4]3 shows the equal-time correlation func- 
tion C c f )C f ) (x^t) = {4>{x, t) 0(0, t)) — (0(0, t)} 2 evaluated with 
TWA. For t/f > 1, it saturates, indicating that the field 
enters a statistical steady state. One can find a decent 
approximation for at large times t ^> f , assuming 
that (j)(x, t) can be represented as a sum of independent 
QBs, at an average density 1/Af in the (x,t)-plane: 



~ x/2 dx f f/2 dt „ 



A/2 



A 



f/2 f 



x(x',x ,to)x(0',XQ,to), (7) 



which describes the core part of fairly well (cf. 
Fig. [4]). Here, a single QB from the 'multibreather' so- 
lution centered at xo enters: x( x: > x o,to) = ®(f — \x — 
x o\)x( x ~ x o,to, k) (6 denotes the heavy-side step func- 
tion). 

This correlation function is directly accessible in ex- 
periments and should distinctively reveal the presence of 
QBs. Moreover, one could simply perform a direct sta- 
tistical analysis of <fi (cf. Fig.[4^i). 

Summary. - We predict the formation of localized 
modulations in the relative phase field, after suddenly 
switching on the tunnel-coupling between a pair of quasi- 
1D condensates. These 'Quasibreathers' grow out of ini- 
tial quantum fluctuations, mimicking processes that are 
crucially important in other areas like cosmology. They 
can be well described by exact analytical solutions of the 
sine-Gordon model. We derived their mean lifetime and 
density after the system reaches a statistical steady state. 
These predicitions are consistent with our numerical sim- 
ulations. An experimental realization, even with present 
setups seems to be within reach. 
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Supplementary Material 

Derivation of Quasibreather Solution - Here, we sketch 
the crucial steps involved in the derivation of the quasi- 
breather solutions discussed in the main text. The details 
of the rather lengthy calculation will be shown elsewhere. 

In general, the Backlund transform allows for the con- 
struction of arbitrary multisoliton waves for the classical 
sine-Gordon equation (SGE) <fi tt — <j) xx +sm(f) = (setting 
for the moment /3 = m = 1). The idea is that starting 
from an arbitrary solution ip of the SGE, one can easily 
show that the field fulfilling the Backlund equations 
(BE's) 



ld_ 

2~dz 



asm 



+ ?A 



2 or a 



(8) 



(9) 



where a G C is an arbitary constant, z = and 
r = ^y^, is a solution as well. Furthermore, it turns out 
that obtained in this manner is just the original solution 
plus an additional soliton (see for instance [T5]). For real 

values of a, the velocity of the added soliton is given by 

2 — i 

v = ^ 2+ j . In the considered case, the known solution ip is 
the homogeneous part of the phase-field, <&(t), perform- 
ing Josephson oscillations according to $ = — sin<l>. As 
the numerical analysis with the TWA shows, slight inho- 
mogeneties in the field are amplified parametrically, even- 
tually leading to the formation of (standing) breather-like 
excitations, which have a finite lifetime. We construct 
exact solutions of the SGE describing standing breathers 
on top of this oscillating 'background'. As it was shown 
in the main publication, it turns out that the solutions 
constructed here are well suited to describe the observed 
localized patterns in the field 0. 

In the absence of the background, ordinary breathers 
of the SGE can be obtained by adding a soliton and an 
antisoliton to the 'vaccum'. Here, we follow closely this 
strategy. For the purpose of creating a single soliton, the 
parameter a can be choosen to be real. However, in or- 
der to create finally a breather out of a soliton-antisoliton 
pair, we must allow for a G C. Being interested in stand- 
ing quasibreathers (breather + background) , it turns out 
that in a first step, formally, one has to add a single 
'soliton' with \a\ = 1 on top of the background. Solv- 
ing Eqs. ( |8|9[ ) with -0 = one obtains a set of solutions 
(j) a (x,i). Fortunately, for the second step (adding the 
corresponding anti-soliton 0^), one can use the fact [15J 
that 



x(x, t) = 4arctan 



a + 

a — j3 



tan 



*(*Xio) 



is a solution of the SGE as well, as long as both a and 
0/3 were obtained from Eqs. ( 8|9 ) with i\) = <I>. The first 



step (finding the solution for a 'soliton' on top of an os- 
cillating background) involves solving the set of coupled 
differential equations Eqs. (8|9). The second step mainly 
consists in matching a and f3 properly (it turns out that 

a = P*). 

Finally, one finds that the quasibreather solutions x 
have the form 



4 / 2 InLF[flj,t;^ 0j A ,<p] 

X = 77 arctan 

p \tan(p 



i + m 2 



$.(ii) 



The argument of the arctan-function corresponds to the 
function Q introduced in the main text. In fact, there are 
two sets of solutions: 



T_ 

f 



tanh[ 



B 



, cos (p > sm 
, cos ip < sin 



(12) 



Whereas the well known, standing SGM-breather is char- 
acterized by a single parameter (its amplitude), the solu- 
tion here depends on three parameters. The amplitude of 
the underlying background oscillations $0, the maximal 
amplitude of the breather (determined by a parameter 
(f G [0,7r/2]), and the relative phase between breather 
and background oscillations, set by Aq. 



The function A(t) in Eq. (12) grows monotonically 
during the temporal evolution of x- It is given by the 
integral equation 



A(t) 



mfsm(2cp) /• t ^,sin 2 [/3$(i / )/2] 



/' 

^0 



\B\Ht>) 



A . (13) 



Furthermore, B(t) = j3m x $ + 2i sin ip sin( / 5$/2) and 
A(t) = 2/- 1 cos(^cos(/3$(t)/2), where / = 2|cos 2 y?- 
sin 2 (/3$ /2)| 1 / 2 . The solutions for cos^ > sin(^ Q /2) 
are periodic in time [with C(t) = tanh _1 (V^4. 2 — 1 — ^4.)]. 
These are termed 'single quasibreathers' and denoted by 
Xif - For costp < sin(y0$o/2), X 1S periodic in space. These 
'multibreather' solutions Xk are labeled by the wavenum- 
ber k = mf/2 [in this case C(t) = tan~ 1 (vC^TT + A)]. 

From now on, we consider f3$o < 1. Averaging dA/dt 
over one Josephson period, one obtains 2 A ~ (5^t + 2Ao, 
with a frequency difference (see main text) 



sn((p,$ ) 



mf cos cp 
2(1 + sin if)' 



(14) 



The 'multibreather' solution Xk has remarkable proper- 
ties (see Fig. [5|. For t —> —00 [A — » —00] it reduces 
to a phonon Sxk with wavenumber fc, superimposed on 
the oscillating background. It is amplified parametrically, 
yielding 



Sxk oc e k ((t) coskx, 



(15) 



where ((t) is a periodic function. The amplification rate 
Y k = 6Cl(ip(k),<&o). As soon as A ^ —1, the charac- 
teristic breather peaks form and exist during a time set 
by Tk = 4/IV The maximum is reached for A ~ 0, 
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Figure 5. a) 'Multibreather' solution Xk with k = mf3^o/^/8 
and f3&o = 0.257T, evaluated at the center of the excitation 
at x = 0. We plot 5\k — Xk — L -1 dxxk (L denotes the 
system size). The inhomogenity 5(f) is amplified parametri- 
cally. The blue, solid lines show the analytical result, b) 
Cut in x-direction for the times indicated in a). One ob- 
serves the transition from initial, sinusoidal inhomogenities 
into the characteristic breather-peaks at large times. The an- 
alytical solution Xk is checked against a direct simulation of 
the sine-Gordon equation with a 4th-order Runge-Kutta al- 
gorithm (with appropriate initial conditions; red, dashed line) 



where one observes a set of QBs with distance 27r/fc. For 
A — » oo, these decay again. In constrast, Xy describes 
a single QB whose amplitude is modulated quasiperiodi- 
cally with period 2ir/5Q(cp, $o). 

Correlation Function. -At long times t/f ^> 1, as 
reported in the main text, one observes a steady state 
characterized by the perpetual creation and annihila- 
tion of breathers. In order to confirm the presence of 
quasibreathers in the experiment, we propose to eval- 
uate the phase-phase correlation function C^{x) = 
(4>(x,t)4)(0,t)) - (cf)(x,t))((f)(0,t)). Note that as we re- 
strict to initial detunings /3$o < 1, we do not observe the 
creation of solitons after the quench and the correspond- 
ing phase slips in (j). Therefore, we do not have to switch 
to correlation functions like (cos /3c/)(x, t) cos j3(j)(x, 0)). 

One can find a decent analytical approximation for 
based on the 'multibreather' solutions discussed above. 
The 'multibreather' solution is strictly periodic in space. 
However, note that this exact solution of the classical 
sine- Gordon is only valid when the breather amplitude 



is large enough compared to the noisy background. We 
assume that each of these quasibreathers centered at X{ 
(their phases are encoded in the time U) is described by 

X(x;x u t z ) = 0(^ -\x- Xi\)e(^ - \ti\)Xk(x - Xi,1®6) 

Here, Xk denotes the 'multibreather' solution with k = k 
and O is the heavy-side step function. For simplicity, we 
pick the special set of solutions with Aq = (at small 
values of (3&o, Xk depends only slightly on the precise rel- 
ative phase between breather and background, Ao, any- 
way). With this choice of Aq, the 'multibreather' solution 
Xk(x,t) reaches its maximal value at t — and x = nX 
with n G Z. After some time, a 'real' quasibreather dives 
into the noisy background and the 'multibreather' solu- 
tion ceases to be valid. While it is hard to determine 
the precise time when this happens, in a reasonable ap- 
proximation, we cut-off x(x]Xi,ti) for \ti\ > f/2. Here, f 
is the 'lifetime' of the 'multibreather' solution (see main 
text). 

In a given experimental run a, (j) a (x,i) ~ 
Sj=i x( x i x ? it? ) w ^h n a denoting the number of 
breathers at time t (neglecting short-wavelength radia- 
tion). It is understood that the distance between neigh- 
boring breathers \xi — — A. Under the assumption 
of independent quasibreathers and with (</>) ^ as well 
as Cfff^Xjt) ~ C ( / )C / ) (x) for t/f ^> 1, one obtains 



1 1, a 

C^{x) ~ £x(z;a?,t?)x(0;a;? ,tf). (17) 

i=l 

The bar denotes the average over experimental runs. 
Note that only a single quasibreather with \xf \ < A/2 
contributes to the average and that the correlation 
function vanishes in this approximation of independent 
breathers for \x\ > A. Therefore, we can replace the av- 
erage over runs by an average over xq = xf and to = tf 
yielding 

C^x) 

I M/2 rr/2 

-it / dx dtx(x;x ,to)x(fy%o,to), (18) 

r ^J-X/2 J -f/2 

which can be evaluated numerically. 



